function f = curve_point(x) 
global pol_cof
global pol_der
global cur_x
global cur_y
p1 = polyval(pol_cof, x);
p2=polyval(pol_der,x);
f=p2*(p1-cur_y)+(x-cur_x);
end